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In multi-band systems, electrons from different orbitals coexist at the Fermi surface. An attrac- 
tive interaction among these quasi-particles gives rise to inter-band or hybrid pairs which eventually 
condense in a superconducting state. These quasi-particles have a natural mismatch of their Fermi 
wave- vectors, SkF, which depends on the strength of the hybridization between their orbitals. The 
existence of this natural scale suggests the possibility of inhomogeneous superconducting ground 
states in these systems, even in the absence of an applied magnetic field. Furthermore, since hy- 
bridization V depends on pressure, this provides an external parameter to control the wave- vectors 
mismatch at the Fermi surface. In this work, we study the phase diagram of a two-dimensional, two- 
band metal with inter-band pairing. We show that as the mismatch between the Fermi wave-vectors 
of the two hybrid bands is reduced, the system presents a normal-to-inhomogeneous superconductor 
quantum phase transition at a critical value of the hybridization V c = Ao. The superconducting 
ground state for V < V c is characterized by a wave- vector with magnitude |q c | = q c = 2Ao/v/. Here 
Ao is the superconducting gap in the homogeneous state and v/ the average Fermi velocity. We 
discuss the nature of the quantum critical point (QCP) at V c and obtain the associated quantum 
critical exponents. 



I. INTRODUCTION 



In a metal, an external magnetic field polarizes the 
quasi-particles giving rise to a mismatch of the Fermi 
wave-vectors of the different spin bands. In the pres- 
ence of an attractive interaction between electrons with 
opposite spins, the ground state of the system is a ho- 
mogeneous superconductor at least for sufficiently small 



mismatch Skp 



kp. In the other limit of very 



large mismatches, we expect the system to be normal 
(non-superconducting) for physically reasonable attrac- 
tive interactions. An important and fruitful question in 
the study of superconductivity is the nature of the zero 
temperature phase diagram as the mismatch is reduced 
in the normal phase or, alternatively, how the BCS phase 
disappears as 5kp increases. A theoretical answer to this 
problem was given by Fulde and Ferrel [l| and Larkin 
and Ovchinikov These authors have shown the ap- 
pearance of an intervening inhomogeneous superconduct- 
ing phase between the homogeneous BCS and the normal 
phases as 5k p increases. This phase generally known as a 
FFLO phase is characterized by a wave- vector dependent 
order parameter which oscillates in space 0- Different 
types of solutions associated with a single wave-vector 
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can be obtained with different ground state energies [J- 
:_3j. Here we consider the original Fulde- Ferrell (FF) or 
helicoidal solutions such that, the order parameter 
is given by A = Aoe iq x . Experimentally, FFLO phases 
have been elusive and there is no conclusive evidence for 
their existence. 

In cold atom systems, two fermionic species with differ- 
ent densities and attractive interactions provides a con- 
trollable model system to search for FFLO phases 
The mismatch can be tuned by varying the number of 
atoms of the different species. These are clean systems 
in which the strength of the interaction between the dif- 
ferent species can also be tuned. In spite of intensive 
search only for nearly one-dimensional systems evidence 
was found for FFLO type of correlations 

In transition metals, rare-earths and actinides metals, 
electrons arising from different orbitals (s,p,d or f) coex- 
ist at the Fermi surface Q • These electrons hybridize and 
the new quasi-particles have different Fermi wave- vectors 
which among other things depend on the strength of 
the hybridization f§4l0l|. These systems and their inter- 
metallic compounds in the case the dominant attractive 
interaction is between quasi-particles in different bands 
are natural candidates for a FFLO ground state, even in 
the absence of a magnetic field. For a two-band metal, 
the mismatch k F — k b F is directly related to hybridization 
[1, ■ Since this can be controlled by pressure, we have 
an external parameter which allows to probe the phase 
diagram and search for FFLO phases. Inter-band pairing 
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may be particularly important in heavy fermion systems 
Q where the dominant Kondo interaction between con- 
duction electrons and the rather local /-electrons pro- 
vides also a mechanism for inter- band attraction [111 ]. 
Notice that FFLO phases induced by pressure avoid all 
the complications associated with orbital effects due to 
an external magnetic field. 

In this paper we study the phase diagram of a two- 
dimensional (2d), two-band system as a function of the 
mismatch between their Fermi wave- vectors controlled by 
the strength of their hybridization. This is relevant for 
nearly two-dimensional materials with cylindrical Fermi 
surfaces. Some of our results are in common with previ- 
ous ones obtained for field induced FFLO phases in 2d 
materials at finite temperatures [HI, [l3|, Ha - tlTj . As in 
these studies we assume the existence of weak inter pla- 
nar interactions so that the mean field BCS approach 
is justified at zero temperature (l3j . However, these 
works ignore the dynamical aspects inherent to the quan- 
tum critical phenomena. We give here a full treatment 
of the quantum phase transition from the normal-to- 
inhomogeneous superconductor as the Fermi wave- vector 
mismatch is reduced from the normal phase. Since this 
T = transition is continuous or second order, it is as- 
sociated with a quantum critical point (QCP) at a crit- 
ical value V c of the hybridization. The instability of the 
normal state that we consider is that for a FF super- 
conducting state characterized by a single wave-vector 
q. This is the first zero temperature instability as V is 
reduced [TH, [TJ] . Our approach yields the full thermo- 
dynamic behavior near the QCP allowing to identify the 
FF or helicoidal superconducting phase from its precur- 
sors effects. 



II. MODEL HAMILTONIAN 

The Hamiltonian of the two-band system with hy 
bridization and an inter-band interaction is given by |8| 



H = H - Mq?1 q + H int 

fc,a 

ka ka 

Hint = 



(1) 



ka 



9 a l<T b -k>-a b -k*a>k-<T, 

k.k' a 



where at CT a^cr {bka) are creation and annihilation 

operators for the a (b) particles and e£ are their disper- 
sion relations, defined by e£ = t% — /Lt a , where = 
and [i a is the chemical potential of the (non-interacting) 
a-particle (a = a, b). Since the inter- band interaction 
between particles a and b is attractive, we take g < 0. 
In the case of a heavy fermion system, the a-band is a 
wide band of conduction electrons and b a narrow band of 



/-electrons with effective mass mb ^> m a . Multi-band su- 
perconductivity was previously studied by Suhl et al. [l8[ 
and Moskalcnko ■ These authors differently from the 
present study considered intra-band interactions, with 
an inter-band term that transfers Cooper pairs between 
the bands. In our case superconductivity is associated 
with the existence of a finite anomalous correlation func- 
tion, < a\ (J b\_ (J >. Since the conduction electrons are 
described by plane waves, which contain all the harmon- 
ics, this correlation function vanishes only in the normal 
phase. 

Next we calculate the zero temperature response of 
the two-band system to a wave- vector and frequency de- 
pendent fictitious field that couples to the superconduct- 
ing order parameter [l(J of interest, < a[ (J b\_ a >. The 
Hamiltonian associated with this coupling is given by 

m, 



-a), (2) 



where the frequency u has a small positive imaginary 
part to guarantee the adiabatic switching on of the ficti- 
tious field. The linear response of the system to this field, 
is obtained using perturbation theory for the retarded 
and advanced Green's functions We consider that 

in the two-band system the hybridization is sufficiently 
large, such that, in this initial condition the system is in 
the normal phase. In this case the superconducting order 
parameter is zero in the absence of the fictitious field h q . 

The superconducting response to the time and in- 
dependent fictitious field, within the BCS decoupling is 
obtained in the form [Toll. 



SA?=x(q,e>)h q 



l-5Xv(9,w) 



(3) 



where, SA q b 



Efc < a k + q «bl k _ a >, and the non 
interacting pair susceptibility Xv(Qi u ) is given by, 

W(4-,) -/(<;*+,) 



e l-q + e t+q ~ W 



(4) 



where 1 and 2 refer to the new hybridized bands (see 

Fig. ED, 

4' 2 = vf{k-k F )i i v, 

where Vf = (vj + Vj)/2 is the average Fermi velocity and 
we took, [i a = kp/2m a . For simplicity we considered a 
situation of perfect nesting for the unhybridized system, 
such that, the two bands cross the Fermi surface at the 
same Uf for V — as shown in Fig. [TJ The mismatch 
of the new Fermi wave- vectors of the hybridized bands is 
given by, SkF = hp — kp = 2V/vf. In the equation for 
the susceptibility we can neglect spin indexes since the 
normal system is paramagnetic. 

In the normal system, at zero temperature, as the 
strength of the hybridization V decreases, the condi- 
tion 1 — gXvili 0) = in Eq. [3] is eventually satisfied. 
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N J »J 




5k F = k F -k 2 F = — 

f r v f 




FIG. 1: (Color online) Simplified band-structure used in the 
calculation of the dynamic susceptibility. For the original 
bands unhybridized bands (dashed lines) we assume perfect 
nesting, such that, they cross the Fermi surface at the same 
wave-vector. The hybridized bands bands are shown as full 
lines and cross the Fermi surface at kp and kp, with a mis- 
match 8k f = kp — kp. 



In this case the interacting pair susceptibility x(q 7 0) 
diverges implying that even in the absence of the fic- 
titious field (h q = 0) there is a spontaneous appear- 
ance of superconductivity in the system. The condition 
1 — g\\? (q, 0) = is a generalized Thouless criterion [2(| 
for a g-dcpcndent superconducting instability. This in- 
stability can be triggered, as discussed below, either by 
decreasing the mismatch between the Fermi surfaces, i.e. 
reducing the hybridization, or by increasing the strength 
of the attractive interaction. The quantum phase tran- 
sition from the normal to the FF or helicoidal super- 
conducting state at the critical value of hybridization V c 
(or g c ) is continuous and in the three dimensional case 
it's universality class has been obtained [l(J. The value 
of the wave- vector q c characterizing the helicoidal super- 
conducting state for V < V c (q) is that for which the 
condition 1 — gXv(q>fy = is first realized. It corre- 
sponds to that where V c (q) determined by the condition 
above has a maximum (V c (q c ) is a maximum). Due to 
the cylindrical symmetry of the problem any linear com- 
bination, A(r) = J2 m A m e !qm r with |q m | = q c has the 
same critical value of hybridization [13] . This degener- 
acy for V < V c is removed by additional interactions or 
non- linear terms in the gap equation (l3l |. Any of these 
phases will have the same quantum critical behavior ob- 
tained below. 



A. The nature of the transition 

The non-interacting particle-particle, or pair dynamic 
susceptibility given by Eq. 0] can be easily obtained at 
zero temperature. It is given by 



Xv («>< 



pin 



2lo c 
V 



(5) 
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FIG. 2: (Color online) The logarithm of the function in Eq.[7] 
This has a maximum for q = 1 with discontinuous derivatives 
at this point. 



where p is the density of states at the Fermi level of the 
2d system, q = Vfq/2V and w c an energy cut-off. The 
real part of the static susceptibility is given by, 



Rext?(g,w = 0) = p 



2uj r 



l n ^£_ m [l + v /7— 
h 



l 2uJc l - 
In — In q 

h 



,q > 1. 



,9<1 
(6) 



The Thouless criterion for the appearance of supercon- 
ductivity is given by, 1 — g$texv c {q, 0) = and coincides 
with the condition for the divergence of the static inter- 
acting pair susceptibility, w = 0) in Eq. [3] 

At zero temperature, for a fixed value of the interaction 
and momentum, as the hybridization V is reduced, the 
normal system becomes unstable to a superconducting 
ground state at a critical value V = V c determined by the 
above condition. Using that the superconducting gap for 
V = can be written as, Ao = 2u} c e~ 1 / gp , we obtain for 
Vc(q), 



V c (q) = 



Ao 



1 + V 7 ! 
A 



;,q< 1 



,q> 1. 



(7) 



The logarithm of this function is plotted in Fig. [21 It is 
continuous and has a a sharp maximum at q = 1. How- 
ever, it's derivative is discontinuous and different from 
zero at the maximum. 

The value of q for which the instability first occurs is 
that for which V c (q) has a maximum, namely, q = q c = 1. 
At this value of q, V c = Aq. Also, using that kp — kp = 
2V/vf, we get q c = kp — kp, i.e., the wave vector of 
the instability is exactly that connecting the two Fermi 
surfaces at kp and kp, as shown in Fig. [3[ This is also 
true in Id, but not in 3d (Toj . 

Next we expand the denominator of the interacting 
dynamic pair susceptibility, Eq. [31 near V = V c , q = q c , 
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FIG. 3: (Color online) The wave-vector q c of the helicoidal 
superconductor is that joining the two Fermi surfaces at kp 
and kp, i.e., q c — k F — kp. 



and small frequencies w. For q/q c > 1 and 1 — 2uj /v/q c < 
q/q c < 1, we get: 




W (8) 



1 -gxv c (q,u) « gp 



This allows us to write the Gaussian part of the action de- 
scribing the normal-helicoidal superconductor quantum 
phase transition [22, HH . It is given by, 



dq / du! 



V-V c 




|A,(?,o;)| s 



(9) 

The quantity S = (V — V c )/V c measures the distance 
to the quantum critical point (QCP) associated with 
the normal-to helicoidal superconductor instability at 
V c = Aq and with q c = 2Aq/i>/. Since hybridization 
increases with pressure (P), near the QCP we can write 
5 = (P — P c ) / P c where P c is the critical pressure below 
which the g-supcrconductor is stable. Within this Gaus- 
sian approximation, we can define a reduced correlation 
length £/a = 1/5, where a is the lattice spacing, such 
that, the zero temperature correlation length exponent 



takes the value v = 1 for this model. Also, we can imme- 
diately identify [22j that the dynamic critical exponent 
z that scales the energy, or frequency, at the QCP is 
given by z = 2. Even at the Gaussian level this quantum 
phase transition is in a new universality class since the 
critical exponents are non-standard Gaussian exponents. 
For example, we find v = 1 instead of the usual Gaussian 
correlation length exponent vq = 1/2. 

The free energy associated with the Gaussian action in 
two dimensions, Eq. [9j is given by [22j . 



2 1 f q ° , 
-JTT^T / dqq 



o e*-l 



tan 




(10) 



It can be written in the scaling form, 

focW^F^], (11) 
with d = 2, v = 1 and z = 2. The characteristic tem- 
perature T* = V C C Z = V c \5\ uz = V C \S\ 2 . In the non- 
critical side of the phase diagram in Fig. 0] there is a 
crossover temperature between two different regimes. For 
T << T* , the free energy is given by, 



/ = 



-C(3/2)T 3 / g 



47r 5/2 



foc-r^U + acO]- (12) 



At low temperatures, T < T*, the specific heat C/T = 
—d 2 f/dT 2 has a power law behavior C oc y/T with 
a coefficient that increases as the QCP is approached 
(£ _1 cx 8 — > 0). So at low temperatures before the insta- 
bility, the normal system behaves as a non-Fermi liquid. 
For T > T*, the specific heat beside the C oc \/T term 
has contribution C/T cx InT. The latter is the scaling 
contribution C/T cx T^ d - Z ^ z which for d = z gives rise 
to a logarithmic term. 

Now let us assume that the action close to the QCP 
has an expansion in powers of the order parameter, 



/ dqj du[xa\q^)] \A(q,u)\ 2 + 

~[dcfi / ]~[ duiu({qi}, {u)i})A(qj, ui)A(q^, q; 2 )A(g3, ^3)A(gl, w^Sj'y] <&)5(y] Uj), (13) 



r 



with i = 1 to 4 and a quartic coefficient u({qi}, {woi})- 
If we neglect all q and u dependences of this term, 
such that, u({q"i} , {u>oi}) = u , simple power counting 
shows that the correlation function of the order parame- 
ter scales as A' 2 = b d+z ~ 2+71 A 2 and we find the critical 
exponent 77 = 1. Furthermore the quartic interaction 



among the fluctuations scales as: 

u< = b 2 -^u , 



(14) 



such that these interactions are irrelevant, in the renor- 
malization group sense, for d+z > 2. Then in the present 
case this implies that the Gaussian action gives the cor- 
rect exponents describing the quantum critical point. 
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However [21[ , uo is dangerously irrelevant for d + z > 2 
implying departures from naive scaling [22j|. For exam- 
ple, the shape of the temperature dependent critical line 
shown in Fig. H is given by, T C (V) oc (V c - V)^ with the 
shift exponent, ip = vz/(l + v9 u ), where 9 U = d + z — 2 is 
the scaling exponent of the dangerously irrelevant quar- 
tic interaction uo [24]]. We obtain tp = z/(d+z — l) = 2/3 
different from the mean-field result ip = 1/2 and from the 
naive scaling prediction ip = vz = 2 [22| . The correlation 
length along the quantum critical trajectory (QCT) in 
Fig. [1 diverges as T — > as £ cx {u T x ^)- v = u" 1 ^ 3 / 2 
j25| , which makes manifest the dangerous irrelevant char- 
acter of Mo . This is different from the naive scaling result 
£ cx T~ x / Z which does not take into account the danger- 
ous irrelevant nature of uq- 

We have given above a complete description of 
the normal-to-inhomogeneous superconductor quantum 
phase transition at the Gaussian level. As we have shown 
this is the correct theory of this quantum phase transi- 
tion for a constant uq quartic interaction. However, the 
assumption of taking u({qi}, {u>oi}) = uq may not be jus- 
tified, particularly in two dimensions. For the case the 
frequency-dependence of u is discarded, but its momen- 
tum dependence is considered, it has been shown that the 
quartic interaction obtained from an expansion of the gap 
equation has a diverges at the wave-vector q c = 1 of the 
superconducting instability [l3l [l7j . At finite tempera- 
tures this divergence is removed and this also may occur 
for finite frequencies or moderate disorder. A possible 
consequence of this anomalous behavior of u is to change 
the order of the quantum phase transition [l| [I?} ■ 

Finally, it is interesting to point out that the cor- 
relation length exponent of the incipient FFLO phase 
(vfflo — 1) being larger than that of any competing 
spin density wave {vsdw = 1/2) [HI, the former should 
override the latter instability. 



III. DISCUSSION AND CONCLUSION 

In multi-orbital metallic systems electrons from dif- 
ferent bands coexist at the Fermi surface. These elec- 
trons have different Fermi wave-vectors and this differ- 
ence provides a natural scale for inhomogeneous super- 
conductivity in the case of inter-band attraction. The 
Fermi wave-vectors mismatch exists even in the absence 
of an external magnetic field and can be tuned by ex- 
ternal pressure which directly changes the hybridization. 
Among the multi-orbital metals, the heavy fermion ma- 
terials appear as strong candidates to exhibit the phe- 
nomenon investigated here. The dominant interaction in 
these compounds is the Kondo coupling between itiner- 
ant and localized electrons and this naturally leads to 
an attractive inter-band interaction. Many of these sys- 
tems have tetragonal structures with cylindrical Fermi 
surfaces for which the results we have obtained directly 
apply. Also, two-dimensions favors the inhomogeneous 
FFLO type of phases [l2j and in addition the mechanism 
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FIG. 4: (Color online) Phase diagram of the nearly two- 
dimensional metal near the normal-to-inhomogeneous super- 
conductor quantum instability at V c . The critical line T c oc 
(V — V c ) 2,/3 . QCT (dashed line) refers to the quantum critical 
trajectory (see text). The line T* marks the crossover between 
two different non-Fermi liquid behaviors, as discussed in the 
text. The homogeneous superconductor part of the phase di- 
agram for small V and its boundary with the helicoidal or FF 
phase are not shown in the figure since they are not studied 
here. 



studied here avoids the complications due to orbital ef- 
fects since it does not require an external magnetic field. 
Heavy fermions are in general very sensitive to external 
pressure and this allows to probe a large region in their 
phase diagrams. For Ce based heavy fermion systems, in 
many cases, at least two superconducting phases appear 
as they are driven with increasing pressure from anti- 
ferromagnetic, superconducting and normal states [26| . 
The FF or helicoidal superconducting state is a candi- 
date for the first superconducting instability as pressure 
is reduced from the normal state. This transition oc- 
curs at a superconducting quantum critical point (SQCP) 
whose properties and universality class we have obtained. 
The theoretical description of this quantum phase tran- 
sition requires considering time and space as fundamen- 
tal dimensions [22j |. This is implied by the uncertainty 
principle which must remain scale invariant leading to 
the inextricability of these quantities. We have charac- 
terized the thermodynamic behavior near the QCP and 
this should be useful for identifying the helicoidal phase 
from its precursor effects in actual systems. 

The peculiar quantum critical behavior of our model is 
a consequence that q c = I is a singular point of the pair 
susceptibility. In many real systems, we can observe con- 
sequences of the structure of Lindhard functions partic- 
ularly in low dimensional systems, as the Kohn anomaly 
and the Peierls instability in one-dimensional conductors 
[27j . It is reasonable than to expect that our results 
should hold for real materials. 

Finally, we point out that we have not explored the 
phase diagram for small mismatches where a homoge- 
neous BCS phase exists up to a critical mismatch. This 
problem has been treated in Refs. [l2l [l51 - [I7| . 
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Note added in proof. Recently, we became aware of the 
paper [28| where they calculate the susceptibility close to 
the quantum normal-FFLO transition. 
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